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We present a combined theoretical and experimental study of the electronic structure and equation of state 
(EOS) of crystalline PdN2. The compound forms above 58 GPa in the pyrite structure and is metastable down 
to 1 1 GPa. We show that the EOS cannot be accurately described within either the local density or generalized 
gradient approximations. The Heyd-Scuseria-Ernzerhof exchange-correlation functional (HSE06), however, 
provides very good agreement with experimental data. We explain the strong pressure dependence of the Ra- 
man intensities in terms of a similar dependence of the calculated band gap, which closes just below 1 1 GPa. At 
this pressure, the HSE06 functional predicts a first-order isostructural transition accompanied by a pronounced 
elastic instability of the longitudinal- acoustic branches that provides the mechanism for the experimentally ob- 
served decomposition. Using an extensive Wannier function analysis, we show that the structural transformation 
is driven by a phase transition of the electronic structure, which is manifested by a discontinuous change in the 
hybridization between Pd-d and N-p electrons as well as a conversion from single to triple bonded nitrogen 
dimers. We argue for the possible existence of a critical point for the isostructural transition, at which massive 
fluctuations in both the electronic as well as the structural degrees of freedom are expected. 
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I. INTRODUCTION 

Solid compounds of some of the platinum metals (PM) and 
nitrogen have recently been synthesized under high pressure 
and temperature These compounds are of fundamental in- 
terest and possibly technological importance due to their un- 
usual or extreme properties. IrN2, for example, has a bulk 
modulus of 428 GPa and is conceivably superhard. 2 Exper- 
imental and theoretical results^ - - indicate that the crystal 
structures of PM-nitrides depend on the group number of the 
parent metal. Specifically, it has been proposed that the ni- 
trides of Os and Ru crystallize in the marcasite structure^ 
Rh2£ and Ir (Refs. d, 0, and|H) in the marcasite and badde- 
leyite structures, respectively, and Pt (Refs. 0,0, andlil-Hoh 
and Pd , 3 i 8 i n in the pyrite structure. First principles calcula- 
tions have suggested alternative low-enthalpy structures that 
have not yet been experimentally observed^ 2 - - — 

One common denominator of all platinum metal nitrides 
studied to date is that compound formation leads to a large 
increase in the separation between the transition-metal ions 
compared to the pure metal phases, which can lead to local- 
ization of the transition-metal d electrons at low pressures. In 
the case of PtN2, IrN2, and OsN2 the bd electrons are ex- 
tended to the degree that significant hybridization with N 2p 
electrons occurs at all pressures. As a result, the behavior 
of these compounds under pressure is uncomplicated; their 
bonding characteristics as well as optical properties remain 
qualitatively unchanged. They are recoverable to ambient 
conditions after synthesis at high pressure and temperature. In 
stark contrast, the nominal PdN2 compound does not appear 
to be recoverable 3 and, as will be shown in this paper, exhibits 
strongly pressure-dependent Raman activity. First-principles 
calculations suggest that both PtN2 and PdN2 are metastable 
at low pressures^ yet only PtN2 has been recovered to am- 
bient conditions. Enthalpy calculations alone thus cannot ex- 
plain the decomposition of PdN2. Bearing in mind that the 



atomic Ad wave functions of Pd have a shorter range than 
their 5d counterparts in, e.g., PtN2, the observed anomalies 
could result from electron localization. It is well known that 
electron localization is the driving force for highly correlated 
electronic behavior, leading to, e.g., high-T superconductivity 
as well as metal-insulator transitions. This kind of behavior is 
most easily seen in 3d- transition-metal compounds. 



In this paper, we present experimental Raman intensities 
which are strongly pressure dependent. We use electronic- 
structure calculations to show that the observed anomalies of 
the PdN2 compound are related to changes in hybridization 
between Pd-4d and N-2p electrons. We find that the stan- 
dard exchange-correlation (XC) functionals within density- 
functional theory, which in the past have so accurately de- 
scribed the 5d platinum metal compounds, fail to reproduce 
the experimental observations. By resorting to a hybrid XC 
functional that provides a much improved description of lo- 
calization, we show that it is possible to explain the origin of 
the apparent decomposition by an isostructural phase transfor- 
mation, corresponding to the conversion from triply to singly 
bonded nitrogen dimers that is accompanied by a pronounced 
long wavelength phonon instability All We furthermore pro- 
vide arguments for the existence of a critical point associated 
with this transition, where strong electronic as well as ionic 
fluctuations should occur. 



This paper is organized as follows. In Sec. HI we discuss 
the experimental and computational methods used in this pa- 
per. The resulting Raman spectra and calculated electronic 
structure as a function of atomic volume, pressure, and inter- 
nal coordinates are addressed in Sec. Jill The implications of 
our findings are discussed in Sec.HVl 
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II. METHODOLOGY 



A. Experimental details 



and Gaussian smearing with a width of 0.1 eV was used to 
determine the occupation numbers. Atomic positions were re- 
laxed with ionic forces converged to 20 meV/A. 



Palladium nitride was synthesized under high pressure and 
temperature using a laser-heated diamond-anvil cell (DAC). 
Raman measurements were carried out using a Jobin Yvon 
HR460 grating spectrometer and a Roper Scientific liquid N2 
cooled charge coupled device (CCD) detector. Spectra were 
acquired using a grating of 1 800 lines/mm, providing a nomi- 
nal resolution of 1 cm -1 . No attempt was made to correct for 
instrumental nonlinearity as this was assumed to be unimpor- 
tant over the frequency range of interest. A continuous wave 
linearly polarized Ar ion laser operating at 488 nm was used 
to excite and measure the Raman spectra.— The laser was fo- 
cused onto the sample using a Mitutoyo 20 x apochromatic 
near-infrared corrected objective providing a focal spot on the 
order of 1 /im. A continuous wave Nd:YAG laser operating 
at 1064nm that was aligned collinearly with the probe laser 
was used to heat the reactants. Diamond-anvil cells equipped 
with diamonds of culet dimension 300 /im and rhenium gas- 
kets were used to generate the necessary pressure. Pressures 
were determined according to the ruby fluorescence scale of 
Mao. 16 Ruby chips within 20 /im of the probed area of the 
sample were used for this purpose and yielded pressures that 
we estimate to be accurate to within 2 GPa. 

Raman spectra were obtained in two experimental runs. 
For Run 1 nitrogen (Air Products, ultrahigh purity grade) was 
loaded at room temperature under high pressure into the DAC 
together with the Pd metal while for Run 2 it was loaded as 
a cryogenic liquid (99.998%, Praxair). The Pd metal was the 
same in both cases (99.95%, 200 mesh, Alfa Aesar, 12068). 
Pd metal was thermally insulated from one diamond by a thin 
layer of ruby. Before introducing liquid nitrogen into the DAC 
cavity the entire DAC was placed in a dry-nitrogen purged 
glove box and heated at 383 K for 12 h after which the sample 
chamber was sealed. The DAC was then removed from the 
glove box, placed in an open-mouth dewar and rapidly im- 
mersed in liquid nitrogen. The dewar mouth was covered to 
generate a positive N2 gas head pressure to discourage atmo- 
spheric contamination (e.g., by H2O or O2). A few minutes 
later the sample chamber was opened to introduce liquid ni- 
trogen and then closed and pressurized^ 



1. Raman intensities 

In a Raman- scattering experiment, the sample is exposed 
to a monochromatic light source. Photons with angular fre- 
quency uj l are absorbed by the crystal and emitted at a num- 
ber of shifted frequencies uji corresponding to creation of 
phonons. The likelihood for such an event to occur is de- 
scribed by second-order perturbation theory and can be ex- 
pressed in terms of Raman susceptibility tensors^ The mo- 
mentum carried by the photons is assumed to be small, and 
therefore only zone-centered phonons are considered. 

In the context of density-functional perturbation theory, 
the Raman susceptibility tensors can be expressed as a sum 
over third derivatives of the total energy with respect to the 
electric field and phonon displacements, and in the case of 
longitudinal-optical phonons, a contribution arising from the 
generated macroscopic electric field. For details we refer to 
the papers by Gonze, Lee, and Veitheni^i"— The tensors were 
computed within the LDA using the ABINIT software* 2 ^ We 
employed a plane- wave kinetic-energy cutoff of 30 hartree 
and sampled the Brillouin zone with alOxlOxlO Monkhorst- 
Pack grid. Norm-conserving Troullier-Martin pseudopoten- 
tials were employed. 28 For each volume the atomic positions 
were relaxed until the forces were less than 5 meV/A. To en- 
able comparison with the unpolarized experimental spectra, 
we averaged the calculated intensities over the individual po- 
larizations and crystal directions. 



2. Wannier functions 

Wannier functions provide a powerful tool to visualize and 
analyze bonding in extended systems^ A general transforma- 
tion between the Bloch and Wannier representations is given 
by 
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B. Computational details 

All first principles calculations, excluding the Raman inten- 
sities, used the projector augmented wave (PAW) method 18 
as implemented in the Vienna ab initio simulation pack- 
age (VASP)i2 and the Heyd-Scuseria-Ernzerhof XC func- 
tional (HSE06). 20 Comparisons are also made with the local- 
density approximation (LDA) and generalized gradient ap- 
proximation (GGA) as parametrized by Perdew, Burke, and 
Ernzerhof. 21 The Brillouin zone integrations were performed 
using 8x8x8 and 6x6x6 Monkhorst-Pack /c-point grids 22 
for structural relaxations and calculation of Raman frequen- 
cies, respectively. The plane- wave cutoff was set to 520 eV, 



where Urnn is a unitary matrix. By defining a measure of the 
spatial spread functional of the Wannier functions 



n = J2 [( w nW\r 2 \w n (0)) - (w n (0)\r\w n (0)Y 



(2) 



which can be discretized in reciprocal space, Marzari and 
Vanderbil^ were able to devise an algorithm to obtain 
so-called maximally localized Wannier functions (MLWFs). 
This algorithm only requires the overlap between the cell- 
periodic Bloch functions of neighboring /c-points 



M mn = (^mk|^nk+b). 
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FIG. 1. (Color online) Measured Raman spectra obtained after re- 
action between Pd and N2. Spectrum obtained in experimental Run 
1 at ~ 57 GPa. The inset shows in more detail the frequency range 
from 750 to 1200 cm -1 . Peaks are indexed assuming pyrite PdN2. 
The low-frequency peaks are due to nitrogen. 



The MLFWs were constructed by interfacing VASP with the 
Wannier90 program package. 31 To compute we use d 
the earlier preconverged charge densities to generate Bloch 
functions on a 6 x 6 x 6 Monkhorst-Pack grid. For the iterative 
minimization to reach the global minimum (with real- valued 
MLWFs) it is important to provide suitable initial guesses and 
their projections onto the Bloch states. These were calculated 
in real space using the soft part of the PAW wave functions. 
All MLWF images presented in the following were rendered 
by using the POV-RAY 32 software. 



III. RESULTS 

A. Raman spectra 

Upon heating the sample, a reaction occurred above ap- 
proximately 58 GPa. This was indicated by the appearance of 
two new and relatively intense peaks in the Raman spectrum 
(A g and E g in Fig. A third weak feature on the high- 
frequency side of the most intense peak was observed (T g in 
Fig. [T]). We obtained a rough estimate of the synthesis tem- 
perature on the basis of Stokes to anti-Stokes ratios of the new 
peaks at the temperature at which they appeared. For the most 
intense peaks at 800 cm -1 and 1050 cm -1 , we found values 
of 814 K and 890 K, respectively. According to group theoret- 
ical arguments, the Raman active modes of the pyrite structure 
are A g , E g , and T g ^ There are three distinct T g modes, and 
therefore five modes that can in principle be observed. Here, 
we surmise that two T g modes are missing in the experimental 
spectrum due to inherent weakness or proximity to one of the 
more intense modes. This is similar to the case of pyrite-PtN2 
where only four modes were observed and a weak T g mode 
was conjectured to be masked by one of the intense E g or A g 
modes 
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FIG. 2. (Color online) Raman spectra for a range of pressures. Left 
panel: experimental unprocessed data from Run 2. The lowest spec- 
trum is obtained from the recovered sample after increasing pressure 
to 36 GPa from a minimum of 1 1 GPa. Right panel: calculated spec- 
tra (LDA). The spectra are normalized to maximum peak height. 



Several Raman spectra were recorded at room temperature 
as a function of decreasing pressure. Figure [2] shows the ex- 
perimental as well as the calculated Raman spectra of PdN2 
for a representative set of pressures. In general, the agreement 
is very good over the entire range, confirming the dominance 
of the A g and the E g modes in pyrite-PdN2 . 

Below 18 GPa we did not obtain any convincing spectra. 
However, a type of spectrum characteristic of PdN2 reap- 
peared upon increasing the pressure from a minimum value 
(specifically from 11 to 36 GPa, see Fig. Q]).— In an addi- 
tional third experiment, very weak spectra were observed at 
room temperature down to 14 GPa while room temperature x- 
ray diffraction spectra were not observed below 13 GPa. This 
strongly suggests that PdN2 decomposes below this pressure. 

While the PtN2 and IrN2 compounds in our earlier 
experiments^ exhibited Raman spectra whose intensities 
were not noted to be strongly pressure dependent, the rel- 
atively intense room-temperature Raman spectrum of PdN2 
weakens significantly as the pressure is lowered. We at- 
tempted to quantify this positive pressure dependence in PdN2 
by keeping the laser power reaching the diamond-anvil cell 
constant^ The resulting pressure dependence of the normal- 
ized spectral intensities, i.e., the Lorentzian areas for the A g 
and E g modes normalized by acquisition time is displayed in 
Fig. [3] It can be well described by an exponential function. 

This behavior is indicative of a significant increase in the 
electromagnetic screening of the incident laser beam by the 
crystal as the pressure is reduced. In the inset of Fig. [3J we 
show the HSE06 band gap of PdN2 as a function of pressure. 
It clearly narrows as the pressure is reduced, leading to an 
exponential increase in the number of carriers excited across 
the gap, which in turn decreases the volume penetrated by the 
laser light. The theoretical K band gap vanishes at about 
1 1 GPa, leading to metal-like optical response of the mate- 
rial. The inset in Fig. [3] also shows the pressure dependence 
of the GGA band gaps of the two other insulating PM nitrides 
synthesized to date, PtN2 and IrN2 . In contrast to PdN2 Ra- 
man spectra of these compounds have been measured over the 
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FIG. 3. (Color online) Experimental Raman intensities of the A g 
(squares) and E g (circles) modes as a function of pressure. The 
solid lines have been fitted to the filled data points only. The in- 
set shows the calculated band gaps of pyrite-PdN2, pyrite-PtN2, and 
baddelyite-IrN2. The band gaps of PtN 2 and IrN 2 are underestimated 
since they were calculated using a GGA exchange-correlation func- 
tional. 
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FIG. 4. (Color online) Equation of state as obtained from theory 
(HSE06, red solid, PBE, green dashed, and LDA blue dotted) and 
experiment (Ref. 3) (yellow dashed-dotted) state. The inset shows 
the nitrogen bond length as a function of atomic volume and the gray 
areas indicates the two-phase region. 



entire experimental pressure range. This is consistent with 
the calculated band gaps which are always nonzero at non- 
negative pressures. 



B. Equation of state and two-phase region 

Figure [H compares equation-of- state (EOS) data from x- 
ray diffraction measurements 3 as well as LDA, GGA, and 
HSE06 calculations. The HSE06 calculations are in excel- 
lent agreement with the available experimental data which 
end at 13 GPa. Below this pressure the HSE06 calculations 
show a two-phase region in which the EOS is multivalued 



FIG. 5. (Color online) Experimental and theoretical HSE06 Raman 
frequencies as a function of volume. The corresponding pressures as 
obtained from the equation of state are shown at the top of the figure. 
The gray area indicates the two-phase region (see Fig. [4]). 



corresponding to an isostructural transition within the pyrite 
structure. As shown by the inset in Fig. [H the transition 
is accompanied by a dramatic change in the nitrogen dimer 
bond length, which suggests a conversion from single to triple 
bonds as the pressure is decreased. We note that a change 
in bonding characteristics is also obtained by employing the 
LDA or GGA XC functionals (see discussions in Refs. [g and 
[TO) . However, these two functionals do not yield a two-phase 
region. 

In the previous section we have related the disappearance 
of the Raman spectrum of palladium nitride at low pressure 
to the closing of the band gap. As can be seen in the in- 
set of Fig. [3j the band gap as a function of pressure is V 
shaped: it closes with increasing volume but reopens upon 
further expansion. If the electronic structure undergoes an 
abrupt transition in this region, one would expect a number 
of physical properties to exhibit unusual behavior. In fact, the 
HSE06-calculated Raman frequencies change very rapidly as 
seen in the gray area in Fig. [5j In particular, the A g mode, 
which corresponds to the stretch mode of the nitrogen dimers 
in the pyrite lattice, first sharply drops and subsequently rises 
with pressure. Even more dramatic effects are observed in the 
long wavelength-acoustic phonon modes. In particular, the 
bulk modulus shown in Fig. [6] becomes negative in this region 
which provides a simple microscopic explanation for the de- 
composition below 13 GPa observed in experiments. 

It needs to be emphasized that unlike the HSE06 results 
the LDA and GGA data predict a transition that involves no 
latent heat, a positive bulk modulus (Fig. [6]), and no unsta- 
ble phonon modest These XC-functionals cannot explain the 
experimental observations which implies that the isostructural 
phase transition in PdN2 involves subtle electronic effects. In 
this context, it is also important to point out that for pyrite- 
PtN2 we observed a shortening of the bond length at nega- 
tive pressures. Furthermore, the relation between volume and 
pressure is single- valued at all pressures, regardless of the XC 
functional (LDA, GGA, and HSE06). This provides a consis- 
tent and sensible explanation for the experimental observation 
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FIG. 6. (Color online) Bulk modulus as obtained from theory 
(HSE06, red solid, and PBE, green dashed). The gray area indicates 
the two-phase region. 




r rxMrR rxMrR r x m r 



FIG. 7. (Color online) Orbital-projected band structures of PdN2 at 
three representative pressures. The line color indicates the relative 
admixture of Pd-d (red) and N-p (green) states. Note the band cross- 
ing at the R point in the middle panel. 



that PdN2 decomposes under pressure whereas PtN2 does not. 



C. Electronic structure of PdN2 

The band structure of pyrite-PdN2 is strongly dependent on 
volume. This is illustrated in Fig. [7] where we display the 
orbital-projected HSE06 band structure at three representa- 
tive pressures. At elevated pressures (left panel in Fig. [7]) the 
valence- and conduction-band edges are dominated by N p 
and Pd d states, respectively. At these pressures, PdN2 is an 
insulator with a direct band gap at the M point. When the lat- 
tice is expanded, the band gap decreases until the two-phase 
region is reached. Here, the system becomes metallic with 
a small overlap of bands around the R point. At even larger 
volumes, after the transition has taken place, the system again 
becomes insulating with a direct band gap positioned at two 
thirds along T — X (right panel in Fig. [7]). The orbital charac- 
ter of the band edges is now reversed with respect to the high 
pressure case. 

The transition in the electronic structure that accompanies 
the isostructural transformation is also illustrated in Fig. [8] 
which shows the density of d electrons projected onto spheres 
of radius 1 .4 A around Pd atoms as a function of energy and 
volume. The white segments in Fig. [8] approximately depict 
the V-shaped band gap as a function of volume with a mini- 
mum at 9.5 A 3 /atom. 37 At this volume a massive transition in 
the character of the electronic states takes place as the Pd-e g 
states migrate across the band gap. In the following section, 
we will show that this migration is directly related to the con- 
version from single to triple nitrogen bonds alluded to above. 



1. Wannier function analysis 

To further elucidate the nature of the electronic structure, 
we have performed an analysis of the MLWFs as a function of 
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FIG. 8. (Color online) Partial density of states (arbitrary units) for 
Pd-d states as a function of volume and energy. The vertical bars 
indicate the two-phase region. 



volume. In Fig.Oa), we show the total spread of the occupied 
MLWFs (red solid line). The total spread rises from its initial 
value at elevated pressures until the system reaches the two- 
phase region. In the metallic phase, no separation of occupied 
and unoccupied states can be made. We can, however, con- 
struct approximate MLWFs by considering all Bloch states up 
to the Fermi-level as fully occupied. This is justified since 
in this region there is only a small volume in k- space around 
R where valence and conduction bands overlap. The spread 
obtained in this fashion rises sharply to a cusp [red dotted 
line in Fig.^a)]. At sufficiently large volume, the system be- 
comes insulating again and the spread continues to decrease. 
Hence, the MLWF spread reproduces the behavior expected 
for a system undergoing an insulator-metal-insulator transi- 
tion. As will be discussed in more detail below, the cusp in 
the MLFW spread functionals separates the regions of single 
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FIG. 9. (Color online) (a) Spread of occupied MLWFs. The red 
solid line indicates the sum of all occupied MLWFs. The letters in 
the legend correspond to the MLWFs depicted in Fig. [Tol (b) Total 
spread of occupied the MLWFs as a function of volume and nitrogen 
dimer bond length. The dashed line separates the nitrogen single and 
triple bonded regions. 



and triple bonded nitrogen. 

To fully characterize the electronic structure, we display the 
occupied and lowest unoccupied MLWFs in Fig. [TQl at high 
and low pressures. All occupied MLFWs were optimized si- 
multaneously, while the conduction band MLWFs were ana- 
lyzed separately. 



from the N 2 -2p7r* orbitals. 



3. Low-pressure region 

At pressures below the transition, the Pd t2 g orbitals are 
still occupied, and are joined by the Pd e 9 states [Fig. [TQlf) 
and[T3g)] which here constitute the Pd-N bonding. The tran- 
sition to triply bonded nitrogen dimers is now obvious: we 
obtain one occupied N s-like state [Fig.QIJh)] for each nitro- 
gen, three identical MLWFs per dimer rotated 120° about the 
N-N axis fFig. [TOtiYl. Finally, the unoccupied MLWFs can 
be identified as N2-2p7r* anti-bonding states. This combined 
description of the N-N bond is analogous to the standard pic- 
ture of an isolated triply bonded nitrogen dimer but modified 
to reflect the threefold rotation axis along the N-N bond. 



4. Single vs. triple bonded nitrogen 



In Sec. 1MB I we have shown that during the course of 
the isostructural transition, the bond length of the nitrogen 
dimers decreases sharply. The Wannier function analysis dis- 
cussed in the previous paragraphs has established that in the 
high- and low-pressure regions the nitrogen dimers are single 
fFig. [T0lb)l and triple bonded fFig. [Toll)!, respectively. We 
will now explore how the character of the N-N bond depends 
on volume and dimer separation. The total occupied MLWF 
spread shown in Fig. Ob) displays a maximum in the metallic 
region along a straight line in the V — d plane. Inspection of 
the occupied Wannier functions reveals that this line also sep- 
arates the nitrogen single and a triple bonded regions. Based 
on the N-N separation, Chen et alJ± concluded that the ni- 
trogen dimers in PdN 2 undergo a transition from single via 
double to triple bonded with decreasing pressure. Figure [9}), 
however, clearly shows the character of the N-N bond to de- 
pend on both volume and dimer separation. Also, we do not 
observe double bonded nitrogen dimers. 



5. Band anticrossing 



2. High-pressure region 

The occupied MLWFs in the high-pressure region can be 
classified as follows. While the site symmetry of a Pd ion 
in pyrite is strictly S6, for the present purpose it can be ap- 
proximated as cubic. In this case the five Pd d orbitals (not 
counting spin degeneracy) split up into two e 9 and three t<i g 
orbitals, the latter being lower in energy. This is verified by 
the MLWF-analysis as we obtain three occupied t<i g states 
per Pd [Fig.QIJa)] and two unoccupied e^-like states per Pd 
rFig.fTOtcOl. The latter MLWFs also have a significant con- 
tribution from N-p z orbitals. The nitrogen dimers are singly 
bonded via the N2-2sa orbitals in Fig.fTOlb). Palladium and 
nitrogen are bonded via the N-sp hybridized orbitals shown 
in Fig.[T3c). These consist to a large extent of states derived 



Knowledge of the rotation matrices Umn m Eq. (Q]) en- 
ables us to decompose the band structure to show the rela- 
tive contribution of each Wannier center. The results for the 
high- and low-pressure cases are displayed in Fig. [TT] which 
complements the site and angular momentum projected band 
structures (Fig. [7]). In both cases, the MLWFs correspond- 
ing to the N-N bonds are associated with bands far below 
the valence-band edge. At high pressure the valence- and 
conduction-band edges are dominated by Pd e g states and N2 
2pir* states, respectively. The band character is reversed in 
the low-pressure phase without any significant rearrangement 
of the bands. However, as the states associated with the band 
edges belong to the same irreducible representation, band an- 
ticrossings are observed both as a function of reciprocal lattice 
vector and volume. 



(a) A (b) O _ (e) _ (f) _ (g) 

(c) if (d) £ X jk*i (h) (i) (j) < 

•v** ••• ^ if 

FIG. 10. (Color online) (Center) Pyrite-PdN2 conventional cell, [(a)-(j)] Isosurfaces of MLWFs orbitals, (a)-(c), and (e)-(i) correspond to 
occupied MLWFs in the high- and low-pressure regions, respectively, (d) and (j) correspond to the lowest unoccupied MLWFs. The green and 
orange surfaces denote positive and negative values, respectively, (a) £2.? -like states (three per Pd), (b) N-N single bond (one per dimer), (c) 
Pd-N bond (three MLWFs per nitrogen), (d) Pd e^-like states (e) Pd tig -like states (three per Pd), [(f) and (g)] Pd e^-like states (one each per 
Pd), (h) N s derived state (one per nitrogen), (i) N-N triple bond (three per dimer, each rotated 120° about the N-N axis), and (j) N2-2p7r*-like 
states. 



The eigenvector exchange associated with the band anti- 
crossing first explains the closing and re-opening of the band 
gap. Second, it leads directly to the migration of Pd-e g 
and N2-2p7r* states across the gap as a function of volume. 
Third, the appearance of N2-2p7r* states in the valence band 
in the high-pressure phase explains the increase in the nitro- 
gen dimer bond length and the transition from triple to single 
bonds. 



IV. DISCUSSION 

In Sec. lIII Al using Raman spectroscopy it was shown that a 
compound of Pd and N is synthesized in a diamond-anvil cell 
at pressures above 58 GPa. Convincing Raman spectra were 
not observed below 18 GPa at room temperature (Fig. 0. In 
an earlier study^ we observed that the x-ray diffraction sig- 
nal due to the face-centered-cubic (fee) Pd metal sublattice 
disappeared at pressures below 13 GPa. In conjunction, these 
observations suggest that PdN2 synthesized at high pressures 
decomposes if the pressure is released. The very good agree- 
ment between experimental data and Raman spectra calcu- 
lated for pyrite-PdN2 furthermore confirms the previous struc- 
tural assignment.^ 

Unlike any of the other PM nitrides investigated so far, the 
Raman intensities for PdN2 decrease exponentially with pres- 
sure (Fig. [3]). This intriguing behavior can be understood by 
inspection of the calculated band gap, which displays a posi- 
tive approximately linear pressure dependence in the respec- 
tive range closing at 1 1 GPa. The band-gap narrowing with 
decreasing pressure leads to an exponential increase in the free 
carrier concentration and dielectric screening, causing an ex- 
ponential decay of the Raman intensities. The strong and non- 
monotonic volume dependence of the band gap foreshadows 
the unusual properties of this compound near the decomposi- 



tion point. 

As shown in Ref . [§1 PdN2 and PtN2 are qualitatively sim- 
ilar with respect to the pressure dependence of formation en- 
thalpies. Nonetheless, PtN2 has been successfully recovered 
to ambient conditions while PdN2 decomposes already at 1 1- 
13 GPa. In Sec. IIIIB1 we investigated the EOS to resolve 
the origin of this difference. Above the decomposition pres- 
sure, excellent agreement was obtained between experimen- 
tal data and calculations using the HSE06 hybrid XC func- 
tional, which mixes the conventional PBE XC functional with 
25% screened exchange. Since the symmetry constraints im- 
posed in the calculations prevent the decomposition of the 
compound, the EOS can be mapped out below the decompo- 
sition pressure. The calculations reveal a first-order isostruc- 
tural phase transition right in the vicinity of the experimental 
decomposition pressure. In the two-phase region (gray bar 
in Fig. [4]), several properties become extremal: the band gap 
closes and reopens (Fig. [3]) and the Raman frequencies ex- 
hibit a pronounced dip (Fig. O. Most importantly, however, 
the bulk modulus becomes negative (Fig.© which readily ex- 
plains the experimentally observed decomposition as being 
triggered by a long-wavelength acoustic phonon instability. 
It is important to note that conventional XC functionals such 
as LDA or GGA do neither agree well with the experimental 
EOS nor do they yield a first-order phase transition associated 
with an instability. This clearly highlights the importance of 
electronic effects in understanding this transition. 

As the first-order transition is isostructural, the two phases 
could be electronically distinct. In fact, a systematic Wannier 
analysis shows that on the high- and low-pressure side of the 
transition the nitrogen dimers are single and triple bonded, re- 
spectively rFig.[TQlb)l. Along the transition pathway, the total 
spread of occupied Wannier functions increases as the two- 
phase region is approached, diverges in the metallic transition 
region, and subsequently drops again. Also during the transi- 
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FIG. 11. (Color online) Band structures for the occupied states in the high- (left panel) and low-pressure cases (right panel), respectively. The 
line color indicates the projection of the eigenstates onto the MLWFs. 
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FIG. 12. (Color online) Schematic illustration of the relation between 
pressure and volume that gives rise to a first-order transition at low 
temperatures and becomes second order as temperature is increased. 
The inset exemplifies a possible corresponding phase diagram. 



At low temperatures, this transition is first order but as 
schematically shown in Fig.[l2]the hysteresis diminishes with 
increasing temperature which eventually leads to the phase 
separation line terminating in a critical point. The occur- 
rence of a critical point allows one to avoid the first-order 
transition by first increasing temperature, then pressure, and 
subsequently reducing temperature again as indicated by the 
dashed lined in Fig. [121 The EOS shown in Fig. [4] suggests 
that the isostructural transition in PdN2 could also terminate 
in a critical point. While it is beyond the scope of the present 
work to answer this question, such behavior would provide an 
intriguing pathway to recover the PdN2 to ambient conditions. 



tion the characters of the valence- and conduction-band states 
are inverted corresponding to a band anti-crossing with eigen- 
vector exchange (Figs.[7]and[8]). We find that both volume and 
the nitrogen dimer bond length determine the character of the 
N-N bond, rendering assignments based solely on the bond 
length unreliable (Fig. 0. Furthermore, based on the Wan- 
nier analysis the nitrogen dimers are never double bonded as 
speculated earlier. 11 

Finally, a discussion of the character of the isostructural 
transition is in order. A well-known example of an isostruc- 
tural phase transition is the fcc-fcc transition in cerium.— 
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